Lattice model for cold and warm swelling of polymers in water 



OS 
On 
On 

o 

Q 

5 1 

o 



■4— » 

o3 



13 
G 
O 

o 



> 
in 
oo 
o 

(N 

On 
OS 

-i— > 
c3 



i 

c 

o 
o 



X 



Pierpaolo BruscoliniB and Lapo CascttiS 
Istituto Nazionale per la Fisica della Materia (INFM) and Dipartimento di Fisica, Politecnico di Torino, 

Corso Duca degli Abruzzi 24, 1-10129 Torino, Italy 
(February 1, 2008) 

We define a lattice model for the interaction of a polymer with water. We solve the model 
in a suitable approximation. In the case of a non-polar homopolymer, for reasonable values of 
the parameters, the polymer is found in a non-compact conformation at low temperature; as the 
temperature grows, there is a sharp transition towards a compact state, then, at higher temperatures, 
the polymer swells again. This behaviour closely reminds that of proteins, that are unfolded at both 
low and high temperatures. 

PACS number(s): 05.20.-y; 05.40.Fb; 61.25.Hq; 87.10.+e 



The theoretical modelling of the behaviour of polymers 
in aqueous solution is a long-standing problem which has 
received considerable attention also in connection with 
the problem of protein folding. It is commonly believed 
[EJ that the organization of water molecules in quasi- 
ordered structures around non-polar monomers plays a 
fundamental role in stabilizing the folded "native" state 
of a protein, namely, the biologically active state. In sim- 
ple models of protein folding, however, water is usually 
taken into account within the definition of the interac- 
tions between aminoacids, rather than explicitly consid- 
ering water-monomer interactions. The rationale for this 
approach is that water behaves as a bad solvent for 
hydrophobic monomers, and this can be modeled by an 
effective Hamiltonian with attractive monomer-monomer 
interactions. However, an "exact" effective Hamiltonian, 
obtained from a partial trace over water degrees of free- 
dom, should in principle depend also on the temperature 
and on the conformation of the polymer. Hence it is in- 
teresting to study simple models where the role of water 
can be considered explicitly. 

Here we introduce a lattice model where each site is 
occupied either by a monomer or by water molecules. 
Water-monomer interactions depend on the state of wa- 
ter; water-water interactions are dealt with in a coarse- 
grained way. To distinguish between the effect of the 
presence of water and any other, monomer-monomer in- 
teractions are neglected. The model can be thought of 
as a generalization of the one proposed by De Los Rios 
and Caldarelli (DLRC) ||. Our generalization is an at- 
tempt to model more carefully the water degrees of free- 
dom, yet keeping the description as simple as to obtain a 
model which can be approached by analytical techniques. 
Specializing to the case of a non-polar homopolymer, we 
introduce a suitable approximation scheme and evaluate 
analytically the partition function, the specific heat and 
the average number of water-polymer contacts, which is 
a measure of the compactness of the polymer conforma- 
tion. We find that, for reasonable values of the param- 
eters, the polymer behaves as follows: at low tempera- 



tures, it is found in a non-compact conformation. As the 
temperature grows, there is a sharp transition towards a 
compact state, associated with a sharp peak in the spe- 
cific heat of the water-polymer system; then, at higher 
temperatures, the polymer swells again, smoothly. This 
behaviour reminds that of proteins, that are unfolded at 
both low and high temperatures [0. 

Let us now define our model. We consider a d- 
dimensional lattice with N s sites and coordination num- 
ber z. Each lattice site is occupied either by a monomer 
of the polymer or by a water "cluster" . The polymer is 
made of N monomers. Each monomer can be either hy- 
drophobic (H) or polar (P), to which cases we associate 
a = 0,1, respectively. The sequence of H's and P's is 
fixed, so that the c's are not dynamical variables. The 
number of water clusters is N w = N s — N. Each water 
cluster is a system of m water molecules, to account for 
the fact that several water molecules are affected by the 
presence of a monomer. For each cluster, the state-space 
has an ordinary part, characterized by x — anc l made 
up of an infinity of states i? = 0,1,2,..., oo, and a spe- 
cial part, with X = 1, containing only the state This 
state is special in the sense that it more favourably inter- 
acts with hydrophobic monomers, as we will see in the 
following. The choices of a discrete spectrum and of a 
unique special state are made for the sake of simplicity. 

The model Hamiltonian is written as a sum of a non- 
interacting water term plus a water-polymer interaction; 
no energy is associated to the polymer alone: 



it' ^watcr 

The first term has the form: 



N w 



E* 



(1) 



(2) 



where E* > 0. For each label j the state dj — repre- 
sents ice, i.e., a completely ordered state, with two hy- 
drogen bonds per molecule. Each water-filled site, rep- 
resenting a cell of m molecules, in the d = state has 
2 m hydrogen-bonds, a fraction of which is buried in the 
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bulk of the cell, while the rest connects neighboring cells; 
we set the energy of this state equal to zero. The states 
(Xj — 0; -dj = 1,2,...) are excited states obtained by 
adding kinetic energy to the water cluster: the hydrogen 
bonds are weakened and, eventually, broken. According 
to this coarse-grained description, which disregards the 
details of water-water interactions, the degeneracy g($) 
of the energy level •& is given by the number of ways 
in which one can assign $ "quanta" of energy to the m 
molecules, i.e., the number of ordered partitions of t9 ob- 
jects in m classes: 5(7?) = ( m+ ^ 1 )- 

The special state (x = 1; *) models the partial or- 
dering that water molecules arc believed to assume in 
presence of a non-polar residue 0]. It is indeed an ex- 
cited states of energy E*, with a peculiar character: the 
m water molecules are thought of as geometrically or- 
dered in the cell, so that hydrogen-bond weakening and 
breaking affects only the surface molecules in the water- 
filled site. An accurate evaluation of the degeneracy g* 
of this state would require a detailed specification of the 
microscopic geometry of water, which is beyond the scope 
of the present analysis; nevertheless a crude estimate of 
it can be given by noticing that, in the presence of a 

d-l 

water- monomer contact, m e ff ~ zm a water molecules 
shall rearrange in an ordered conformation to maximize 
hydrogen-bonding. Then, reasoning as before, g* will be 
given by the number of ways to assign $* = E* quanta 
to the TO cff molecules, i.e., g* = ( [[mo ™t' 5 * _1 ) , where [[•]] 
is the integer part of •. 

The second term in the Hamiltonian ([!]) is 

JV 

H mt = J2 A V <Ti[JXj+K(l- Xj )} , (3) 
i=l 

where Ay = 1 if i and j are nearest-neighbors and 
Aij — otherwise, K > (resp. J < 0) is the energy 
cost (resp. gain) of a contact between a H monomer and 
a water site in an ordinary state (resp. special *-state). 
The form of Eq. (^) can be understood if we consider 
a "droplet" of hydrophobic monomers in a water-filled 
lattice, and analyze the energy balance in exchanging a 
water site with a monomer (both taken from the respec- 
tive bulks). If the water is in an ordinary state, each 

d— 1 

contact with a monomer will break up to ~ m d hy- 
drogen bonds between neighbouring water cells. If a > 
is the energy involved in a bond (in units of the spacing 
between the levels), this process yields an energy cost 

d-l 

of K ~ m~3~a, while the cost will be zero if the water 
is in the special *-statc. At the same time, non-polar 
residues interact more favourably with water than with 
each other, due to the permanent dipole moment of wa- 
ter Q . This is taken into account by assuming an energy 
gain J (—a < J < 0) for a contact between water in the 
*-state and a hydrophobic monomer. It should be clear 
that this theoretical framework is not meant to provide a 



detailed microscopic description of the physics of water, 
but rather to account for the basic ingredients of non- 
polar solvation: the existence of an ice-like ground state 
of zero entropy for water, and of a set of excited states, 
some of which may be particularly suited for interaction 
with non polar solutes, but involve a substantial entropy 
loss. 

Let us now study the equilibrium thermodynamical 
properties of the model. To compute the partition func- 
tion it is useful to introduce some notations. The maxi- 
mum number of contacts between water and H monomers 
is obtained when the polymer is in an extended config- 
uration, so that M — (z — 2)Nh + &\ + ctjv, where Nh 
is the number of non-polar residues. Given a configura- 
tion C of the polymer, the energy of the system can be 
written as 

n{C) = g [(1 - x 3 ) + K£j) + X j (E* + M 5 )] , (4) 

where £j = £j (C) is the number of contacts between the 
j-th water cluster and H monomers. Denoting by ni(C) 
the number of water clusters that have £ contacts with 
H monomers, the partition function of the model can be 
written as 

Z = Y. = E C(K(C)}; N) Z{{ne(C)}) ; (5a) 

C {n t } 

Z({n e (C)})=J2 e ~ m{C) > ( 5b ) 

where (({ni(C)}; N) is the number of equivalent confor- 
mations for a polymer of length N, i.e., of the conforma- 
tions with the same set of ng's. Now, to ease the notation 
let us put Z{C) = Z({n t (C)}). Using Eq. (@), and the 
fact that water cells do not interact with each other, so 
that it is possible to factor Z(C) according to the labels 
,7 = 1,..., N w , Z(C) can be written as 

Z(C) = P Q N -f[xT iC} , (6a) 

where xg = Pi/Pq, and 

Pi = g*e-^*+^ + e -PK' 1 . ( 6 b) 

(1 — e p ) 

The mean energy and the specific heat are evaluated 
in the standard way as U = (H) = —d\ogZ/d(3 and 
Cy = dll/dT. The average number of water sites pre- 
senting £ contacts with the polymer can be computed 
as (ni) = xg{d\ogZ / dxg), from which one can ob- 
tain the average (total) number of H- water contacts, 
( n c) = 'Ylit£{ n t)i which is a measure of the compact- 
ness of the polymer. Another interesting observable is 
the average number of contacts with water sites in the 
"-state, given by (n*) = -TdlogZ/dJ. 
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Notice that, even in the simple hydrophobic homopoly- 
mer case, our model can not be mapped onto a polymer 
model with monomer-monomer contact interactions, be- 
cause the partition function (|5|) can not be written in the 
form 



2 = ^C(n H H)e" 



/3/iriHH 



(7) 



where tt-hh = (M — ^2 ^ n f) /2 is the number of monomer- 
monomer contacts and ((jihh) is the number of polymer 
configurations with 7Jhh internal contacts, with h a true 
coupling constant, independent of T and {n{\. 

Before going further with the discussion of the thermo- 
dynamical properties of our model, let us note that the 
partition function of the DLRC model || can be obtained 
from the present one with the substitutions g* — ► 1, 
E* -> 0, (1- e-Py m -> q - 1, due to the fact that 
in the DLRC model all the states are equivalent (with 
zero energy) for pure water, and when a water-monomer 
interaction takes place, the special state has zero entropy, 
while the degeneracy of the excited states is independent 
of the temperature. In both models one has to perform a 
sum over the conformations of the polymer to calculate 
the thermodynamical functions. This is a hard numeri- 
cal task when the length of the polymer grows: for this 
reason we introduce an approximation that allows us to 
evaluate them analytically. First, we test our approxi- 
mation against the exact numerical results obtained by 
DLRC in the case of a homopolymer (with N < 25) on 
a 2-d Manhattan lattice; then we apply it to our model. 

From now on we specialize to the case of a non-polar 
homopolymer of length JVona 2-d square lattice {z = 4). 
Let us now come to the main point of our approach, i.e., 
let us introduce an approximation for the unknown quan- 
tity £({n;(C)}; N) in Eq. (||), representing the number 
of polymer conformations characterized by the same set 
{n-e}. We observe that this number must be strongly de- 
pendent on n c , the total number of water-monomer con- 
tacts, that in turn is related to the number of monomer- 
monomer ?ihh contacts by n c = M — 2nnH, and hence is 
a rough measure of how compact the polymer is. Then, 
our approximation is built in the following three steps: 
1. we assume that (({ni(C)}; N) depends essentially on 
n c and not on each of the rii . More precisely, we assume 
that, at fixed number of contacts n c , the number of con- 
formations in which there are water sites with three or 
more contacts with the polymer is negligible (single and 
double contacts are much more likely to appear, due to 
geometrical reasons), so that C,{{ng (C)}; N) does not de- 
pend on 713 and ri.4. Moreover, we assume also that £ is 
the same for each m, n<i such that ri\ + 2n 2 = n c , i.e., 



C({n e (C)};N) = £(n c ;N) = CV(n HH ) , 



(8) 



where C/v(nHH) is the number of walks of length N with 
"hh = (M — n c )/2 internal contacts. 2. We consider 



((n c ; N) as characterized by two regimes, a globular and 
an extended one, referred to as ( c and ( g , respectively, 
which come into play depending on the value of n c : 



t({ni(C)}) 



C c (n c ; N) if n c > n c ; 
C g (n c ; N) if n c < n c . 



(9) 



This is justified by the fact that the polymer shows, at 
different temperatures, either a globular, compact phase 
or an extended one. The latter involves a high number 
of water-monomer contacts, n c , while the former is char- 

d— 1 

acterized by small values of n c , namely, n c oc N a , as 
in the case of a maximally compact state. We assume 
that a step function in the number of contacts separates 
these regimes; the position of the step, h c , is at present 
unknown. 3. We consider the extended conformations 
as self avoiding walks (SAWs) and the compact ones as 
Hamiltonian walks (HWs): their numbers are both ex- 
ponential in N, and are related by 

( C (N) ~ exp(aA%(A0 , (10a) 

where a is given by 

a = log /j S aw - log /Jhw (10b) 

and /is Aw> Mhw are the respective connectivity constants 
§• 

Resorting to the above assumptions we can evaluate 
the partition function. In the homopolymer case n c is 
always even, so that we set k = n c /2, k mm — n™ ln /2, 
where n" 

k = n c /2 and sum over k, obtaining 

Z = P$ r ~Cg(N) ^(xi,x 2 ;fc min ,fc) 

„ceJV 



z v N is the minimum number of contacts, 



+ e alv iP(x 1 ,x 2 ;k + l,M/2) , 



(11a) 



where -0 is given by 



~>P(x,y,p, q) 



,2p 



{x 2 + x 2(q - p+2) {y - 1) - y [1 + {x 2 - l)y q - p+1 ] } . (lib) 



The thermodynamical observables, U, Cy, (ni) and (n*) 
are then computed according to their definitions. 

The approximation we introduced contains the thresh- 
old parameter, n c , which is not a free parameter, but 
is unknown; its evaluation would require the knowledge 
of the number of walks of arbitrary length and number 
of internal contacts on a lattice. In order to circumvent 
this difficulty, we fix h c applying our approximation to 
the DLRC model, for which we know the results from 
exact enumeration on a 2-d Manhattan lattice || . Using 
our method, the best approximation of the DLRC model 
specific heat is obtained using with n c = n™ m 23/20; we 
choose this value for n c . For the parameter a, Eq. (10b) 



3 



holds, where, for a Manhattan lattice, fJ-sAW — 1-7335 
and ^hw ~ exp(G/7r) = 1.3385 (G is the Catalan 
constant) M. 

Let us now move to our model. Since the discretization 
of the water energy levels is artificial, the ratio a of the 
hydrogen bond energy to the level spacing is a free pa- 
rameter. We choose a = 100; yet, we verified that even 
remarkable changes in a do not imply major modifica- 
tions of the behaviour of the thermodynamic quantities. 
The values of the other parameters have been fixed ac- 
cording to their physical meaning, i.e., we have chosen 
to = 20, E* = 2a, K = a^/m. Given the length N 
of the polymer, we set N w = M = 2N + 2 (which is 
the maximum number of sites which can wet the chain). 
If J < —E*/2, at low T the polymer is extended and 
completely wet by the solvent, as should be a neighbor- 
avoiding walk. Then, raising the temperature above a 
certain T = Tf ~ 13, the polymer becomes compact ((n c ) 
drops to that of a maximally compact conformation), and 
finally, as T grows further, it swells again smoothly (see 
Fig. Q). The specific heat has a peak at Tf and another 
one at a higher temperature, when the polymer is com- 
pact. The latter is related to the excitation of the water 
sites around the polymer in compact conformation out 
of the *-state, as witnessed by the drop of (n*) to zero 
in correspondence of the peak; the former peak at Tf, 
whose height grows with the length N (see Fig. ||), could 
be related to a true phase transition, reminding cold un- 
folding in proteins. The thermal swelling of the polymer 
is present here, but does not have the characteristic of a 
phase transition, as it happens, for instance, in |J. Yet 
this is not a surprise, due to the absence of any kind 
of imposed cooperativity in our model. It is interesting 
to notice that, at low temperatures, the polymer is ex- 
panded even for moderate values of J, i.e., J w — E* /2: 
this means that, in the presence of a monomer, water 
prefers to stay in the excited special state, even when 
the energetic gain J does not apparently compensate the 
energy loss E* . For values of J closer to zero the cold 
collapse is lost: at low temperatures the polymer is com- 
pact, then it swells smoothly at high temperatures. 

To summarize, we have defined and discussed a lat- 
tice model for water-polymer interaction, where water is 
explicitly considered, and introduced an approximation 
scheme that allowed us to analytically evaluate the rel- 
evant thermodynamical averages. We remark that our 
model can not be mapped onto a model with monomer- 
monomer contact interactions alone. We addressed the 
homopolymer case and observed that, for reasonable val- 
ues of the parameters, the polymer is in a compact con- 
formation at intermediate temperatures, and swells when 
the temperature is lowered as well as when it is raised. 
This recalls similar proteins' behaviour. While the "cold 
unfolding" process is sharp, and possibly represents a 
phase transition as N — > oo, thermal swelling is smooth. 
The absence of a sharp transition at high temperature 



could possibly be attributed to the lack of specificity in 
interactions: we expect that the heteropolymer case (to 
be studied next) will show more similarities to real pro- 
teins. The fact that our results are indeed qualitatively 
similar to those of DLRC || , in spite of the differences 
in the model studied (many states vs. two, temperature- 
dependent degeneracy, etc.), suggests that the polymer 
behaviour we both find out is not a peculiarity of a par- 
ticular model, but a property of a class of them: it is 
likely that the crucial thing is to take into account ex- 
plicitly the degrees of freedom of the solvent, though in 
a simplified way. Indeed it appears (see e.g. Ref. p0[ , 
where a model of random heteropolymers is introduced 
and studied) that if the solvent's degrees are neglected 
from the beginning, and hydrophobicity is attributed to 
monomers like a "charge", cold unfolding will not be 
found for non-polar homopolymers, and probably even 
for quenched random heteropolymers. 
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FIG. 1. Specific heat c = CV/(mM) (continuous line, 
rescaled by a factor of 0.2 for graphical reasons), average num- 
ber of water- monomer contacts (n c /M) (dash-dotted), and 
average number of contacts with water in the *-state (n*/M) 
(dashed) as a function of T. Here N — 25 and J = —E*/2. 
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FIG. 2. c = Cv/ijnM) for N = 16 (dashed) and N = 100 
(continuous). J = —E*/2. 
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